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SUMMARY 

The cooling requirements of a plug -nozzle system for a supersonic -cruise aircraft 
engine have been studied. The engine selected provides a thrust of 267 000 newtons 
(60 000 lbf) and is designed for supersonic flight at speeds up to Mach 2.7. The plug is 
sting supported from the turbine frame. The conical plug selected has a 10° half angle 
and may be convection cooled on the sting support and plug up to the 50 -percent point. 

The convection cooling air is then injected as a film over a conical plug extention that is 
truncated at the 75 -percent point on the plug (i. e. , 75 percent of the distance between the 
primary nozzle throat and the tip of the plug). The flight profile selected includes max- 
imum afterburning from takeoff to Mach 2.7. A low afterburner setting is assumed for 
the supersonic -cruise portion of the flight. 

The calculations indicate that, for maximum afterburning, about 2 percent of the en- 
gine primary air, removed after the second stage of a nine -stage compressor, should 
adequately cool the plug and sting support. During supersonic -cruise operation, ram air 
should adequately cool the plug. If about 1. 5 percent ram air is used, the maximum plug 
wall temperature should be about 978 K (1760° R). If 0. 5 percent ram air is used, the 
maximum plug wall temperature should be about 1133 K (2040° R). 


INTRODUCTION 

The cooling requirements have been studied for an air-cooled plug-nozzle system on 
a supersonic -cruise afterburning turbojet aircraft engine to determine the source of cool- 
ing air necessary to minimize the cycle penalty. The airplane chosen for this study was 
assumed to have a takeoff gross weight of 340 000 kilograms (750 000 lbm) and a payload 
of 6. 5 percent of the takeoff weight. (The International System of Units will be used 
throughout this report as the primary system; U.S. Customary units will be included in 



parentheses. Calculations were made using the U. S. Customary units). Each of four 
afterburning turbojet engines were assumed to provide a thrust of 267 000 newtons 
(60 000 lbf), and the cruise flight speed was assumed to be Mach 2. 7. 

In reference 1, the performance of a plug-nozzle system was compared with ejector 
nozzle systems. The optimum performance of a plug nozzle (assuming an uncooled plug) 
and of a variable-flap ejector was excellent and about equal. The performance of the 
auxiliary inlet ejector was slightly lower. However, the plug nozzle has certain advan- 
tages over the ejector nozzles: 

(1) Leakage between moveable surfaces would be significantly lower since the length 
of seal required could be reduced from about 91.4 to 9. 14 meters (300 to 30 ft). 

(2) The mechanisms could be simpler and more durable. 

(3) Jet noise tests indicate that the plug nozzle is a little quieter than the ejector 
nozzles (ref. 2). 

Tending to offset these advantages, however, is the fact that the plug nozzle must be 
cooled during afterburning operation since it is immersed in the high-temperature pri- 
mary jet. Reference 1 showed the effect of nozzle cooling on a typical supersonic -cruise 
mission for (1) a plug nozzle cooled with engine fuel and (2) an air-cooled plug nozzle us- 
ing high-pressure compressor discharge air for the coolant. If engine fuel could be used 
to cool the plug, thus not reducing the engine cycle efficiency, the plug nozzle and 
variable -flap ejector would provide about the same range (about a 3 percent greater 
range than an auxiliary -inlet ejector). However, if compressor discharge air were used 
to cool the plug, the cycle efficiency would drop so that the plug showed little advantage 
in range over the auxiliary -inlet ejector. The purpose of this report is to show that cool- 
ing air for the plug may be removed from an earlier stage of the engine compressor 
without reducing the engine cycle efficiency as much. If this can be done, the plug- 
nozzle system will approach the same range as the variable -flap ejector. 

There are, of course, several trade-offs involved in the design of a nozzle system. 
Among these are nozzle thrust performance, nozzle weight, noise, infrared radiation 
suppression, and some life criterion. Nozzle thrust performance and nozzle weight are 
functions of the type of nozzle and the nozzle -cooling requirements. For a nozzle cooled 
partially by film cooling, the location of the film-cooling slot affects nozzle performance, 
weight, and, of course, the nozzle -cooling requirements. Noise becomes important for 
commercial supersonic -cruise aircraft during flight over land and near airports. Infra- 
red radiation suppression is desirable in certain military applications. Finally, nozzle 
life criterion may be different for commercial or military applications. 

Since the plug -nozzle system is quieter than ejector nozzles and can be made to sup- 
press infrared radiation (by overcooling, for example), a study of the plug cooling re- 
quirements and the resulting performance, weight, and life appears warranted. 

The cooling technique used in this study is similar to the method used in reference 3. 
The plug in reference 3 was strut supported to the engine nacelle. A sting -supported 
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plug was selected for this study to eliminate the struts that were immersed in the hottest 
region of the hot gas. The cooling technique combines parallel -flow convection cooling 
and film cooling. The convection-cooling technique was studied in reference 3; thrust 
performance of an uncooled version of the plug used in reference 3 was presented in ref- 
erence 4. Jeracki (ref. 5) studied the effect of discharging the convection cooling air as 
a film on a plug extension. Finally, Chenoweth (ref. 6) has studied the heat -transfer 
characteristics of a film-cooled plug. Each of these reports has been drawn upon in this 
report in an attempt to optimize the design. 

This report studies the cooling requirements of only the plug and sting support. It is 
assumed that the cooling requirements of the secondary shroud, iris primary, and after- 
burner liner are about the same as an ejector nozzle. 


SUPERSONIC-CRUISE AIRCRAFT FLIGHT REGIME 

The flight profile selected for this study is shown in figure 1 and is assumed to be 
typical of a commercial supersonic-cruise mission. Figure 1(a) presents altitude, and 
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Figure 1. - Supersonic cruise aircraft flight profile. 
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TABLE I. - NOZZLE OPERATING CONDITIONS 


Mach 
number , 
M 

Altitude 

Primary flow 
rate, w 

Nozzle total 
pressure, 

P 8- 

atm 

Nozzle total temperature 

Ambient 

pressure, 

P<r 

atm 

Nozzle 

pressure 

ratio, 

VPo 

— 

\ m 

ft 

K 

°R 

kg/sec 

lbm/sec 

0 

r\ 

\J 

0 

287 

633 

3.03 

1945 

3500 

1.0 

3. 03 

.8 

7 010 

23 000 

189 

416 

2. 10 

1945 

3500 

.405 

5. 18 

1.0 

9 450 

31 000 

172 

378 

1.83 

1945 

3500 

.284 

6.46 

1.2 

10 500 

34 500 

170 

375 

1.87 

1945 

3500 

.241 

7. 75 

1.6 

12 650 

41 500 

166 

365 

1.89 

1945 

3500 

.173 

10. 88 

2.0 

14 500 

47 500 

173 

381 

1.94 

1945 

3500 

.129 

14.99 

2.4 

16 600 

54 500 

182 

401 

1.97 

1945 

3500 

.093 

21.28 

a 2. 7 

18 300 

60 000 

186 

410 

1.97 

1945 

3500 

.071 

27. 6 

b 2. 7 

18 300 - 20 750 

60 000 - 68 050 

186 

410 

1.97 - 1.35 

1347 - 1286 

2425 - 2314 

0.071 - 0.049 

27. 6 


a End of acceleration. 

b Start (first value) and end of supersonic cruise. 


figure 1(b) shows nozzle pressure ratio, both as functions of Mach number, during the 
acceleration portion of the assumed flight. A Brequet cruise was also assumed; the alti- 
tude of the aircraft increases from 18 300 to 20 750 meters (60 000 to 68 050 ft) as the 
fuel is burned. Table I presents nozzle operating conditions used in this study. (Symbols 
are defined in appendix A). Maximum afterburning is assumed during acceleration up to 
Mach 2. 7; the afterburner is assumed to be on at a low setting during supersonic cruise. 
The compressor is assumed to have nine stages. Figures 1(c) and (d) show the 
compressor -inlet and compressor-outlet total temperature and total pressure, respec- 
tively, as functions of Mach number. 


NOZZLE CONFIGURATION 

The afterburner -nozzle configuration is shown schematically in figure 2. The 10° 
plug is sting mounted to the turbine frame. The outer afterburner wall is shielded from 
the afterburner flame with a perforated liner. Turbine discharge air is forced between 
the wall and the liner and film cools the liner. An iris primary nozzle was selected to 
provide variable primary area for the various engine conditions. The iris primary is 
film cooled with the remaining turbine discharge air at the end of the afterburner liner. 

A translating, cylindrical secondary shroud was assumed to vary the hot -gas expan- 
sion ratio and maintain efficient nozzle performance over the required range of nozzle 
pressure ratios. Reference 7 showed that little loss in performance resulted by using a 
two -position shroud. A retracted shroud could be used for takeoff and acceleration up to 
Mach 1.2; then the shroud could be fully extended for speeds up to and including Mach 
2. 7. For the nozzle selected for this study, the retracted position is the -14. 4 -percent 
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Translating shroud 



point (upstream of the primary throat), and the extended position is the 38-percent point 
(see fig. 3). 

The 75 -percent truncation point was selected to reduce the overall length of the noz- 
zle; reference 4 indicates that this truncation should not affect the thrust at high pressure 
ratios but would decrease the thrust coefficient about 1/2 percent at takeoff. 

Figure 3 shows the detailed geometry assumed for this study. The dimensions shown 
for the primary nozzle are for the maximum afterburning position. For nonafterburning 
or low afterburning, the primary nozzle area must be decreased. 

The proposed sting and plug may be convectively cooled to the 50-percent point along 
the plug surface as shown in figure 3, that is , half-way between the primary throat and 
the end of a full cone. The conical extension between the 50-percent point and the 75- 
percent point then could be film cooled with the air discharging from the convective- 


Translating shroud- 

retracted -14. 4 percent ^ Extended 38 percent 



Figure 3. - Plug nozzle geometry (maximum afterburning). Convection cooling to 50 percent, film-cooling to 75 percent. 
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cooling passages. The trade-offs involved in selecting the optimum point to begin film 
cooling will be discussed in a later section of this report. 

Details of the convective -cooling passages in the plug and sting support are shown in 
figures 4 and 5. This method was used successfully on an air-cooled plug in a J-85- 
GE-13 engine (ref. 3). The plug used in the experimental study of reference 3 was 40. 6 
centimeters (16 in.) at its maximum diameter and was convectively cooled to the 60- 
percent truncation point. The plug of reference 3 was strut supported. The fabrication 
technique involved brazing nickel fins onto the inside surface of the high strength plug 
wall. The cooling channel was completed by attaching a thin inner wall. The plug was 
built in axial sections. An expansion seal was incorporated between each section at the 
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(b) Supersonic region. 

Figure 4. - Convection cooling passages. 



0.3175 cm (1/8 in.) 







inner wall, and the outer wall of each section was seam welded to complete the plug. The 
expansion seals are required to allow thermal expansion of the hot outer wall relative to 
the cooler inner wall. 

The fin geometry selected for this study was scaled up by a factor of three from the 
plug used in reference 3; the outer wall was only doubled, however. A preliminary 
stress analysis of this wall is discussed in the section OTHER CONSIDERATIONS. 

The spacing of the nickel fins, and hence the number of fins around the circumfer- 
ence of the plug (fig. 4), was chosen to yield approximately uniform plug wall tempera- 
tures. Uniform high wall temperatures would result in the lowest heat transfer from the 
hot gas to the wall, and, hence, require the least cooling air. Figure 5(a) shows a cross 
section through the support tube, and figure 5(b) shows a section at the maximum plug 
diameter. The dimensions shown for wall thickness, fin thickness, and fin height are as- 
sumed to be the same for the plug and sting support. Fins are shown in the support tube 
spaced at 15° increments. The fins in the support tube are included to maintain uniform 
coolant-passage heights in the support tube and to increase the rigidity of the support 
tube. More fins may be added in the support tube, if desired, to lower the wall tempera- 
tures and increase the strength of the sting support. 

In each section of the plug assembly, the inner wall will be attached to the C -shaped 
fins on the downstream end of that section only. The L-shaped fins are not attached to 
the inner wall. The plug inner cavity should be maintained at the coolant -inlet pressure, 
thus insuring that the inner wall will be pressed against the fins. 


HEAT -TRANSFER DESIGN PROCEDURE 
Convective Cooling 

The procedure used to design the convectively cooled plug reported in reference 3 
was also used in this study. Briefly, beginning with known hot -gas conditions, geometry, 
and coolant -inlet conditions, a heat balance is performed at many stations along the sur- 
face of the plug (see fig. 6(a)). Convective heat transfer from the hot gas to the plug wall 
is included, as is radiation from the hot flame to the wall. Radiation interchange between 
the plug and the afterburner liner, primary nozzle flaps, secondary shroud, and atmos- 
phere were neglected. High-speed flow effects were included. Pressure calculations 
were performed from the end of the convection -cooling passage upstream to the inlet to 
the sting support. Momentum and friction pressure drop were included. The effects of 
these assumptions and equations used are discussed in detail in appendix B. 
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(a) Convection (q c e q conv + q r> f] and q - 0 at end of fin). 


hot gas to wall radiation 


q 2 wall to coolant forced convection 


Primary hot gas 


q 3 wall to ambient radiation 



(b) Film cooling (qj c q? + q 3 *. 
Figure 6. - Heat-transfer model. 

Film Cooling 


The film-cooling design procedure used to predict wall temperatures on the film- 
cooled extension used a modification of the method presented in reference 8, which pre- 
dicted temperatures on a cylindrical ejector. Reference 8 used a wall heat balance which 
included hot -gas -to -wall radiation, convective heat transfer between the coolant stream 
and the wall, internal wall radiation from the ejector exit, external wall-to -ambient 
radiation, and external free convection. The plug wall heat balance used herein includes 
only hot -gas -to -wall radiation, convective heat transfer between the coolant stream and 
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the wall, and external radiation (see fig. 6(b)). 

Flow fields were determined by assuming that the hot -gas jet and the coolant ex- 
panded isentropically along the plug and that the static pressure was uniform across the 
nozzle cross section. The film coolant was assumed to flow between the plug and the hot 
gas. Coolant stream properties beyond the slot were determined by assuming that the 
coolant stream had a recovery temperature equal to the adiabatic wall temperature as 
predicted by the Hatch-Papell (ref. 9) semi -empirical film-cooling correlation, thus ac- 
counting for thermal mixing of the hot -gas and coolant stream. The state of the coolant 
at the slot was assumed to be the state of the fluid at the end of the convective -cooling 
passages. The effect of the fins in the slot on the film-cooling effectiveness was 
neglected. 

The technique for predicting the adiabatic (insulated) wall temperature on a plug noz- 
zle was presented in reference 6 where it was compared with relatively low-temperature 
data. Reference 6 showed that, at takeoff and supersonic -cruise nozzle pressure ratios, 
the correlation provides good agreement between predicted and measured wall 
temperatures. 

The method used to predict plug extension surface temperatures is discussed in de- 
tail in appendix C . 


RESULTS AND DISCUSSION 

Several flight conditions were analyzed in this study. Table I presents the nozzle 
operating conditions assumed for this study. The maximum -afterburning, sea-level- 
takeoff engine condition is critical since, at this condition, no ram air is available for 
cooling the plug. Therefore, the plug cooling air must be obtained from the compressor. 
On the other hand, at the supersonic -cruise flight condition with low afterburning, ram 
air is available and could be used with very little penalty to the cycle efficiency. The 
sea -level -takeoff and supersonic -cruise conditions are examined in detail, and a brief 
look is taken at some of the intermediate maximum -acceleration conditions: Mach 0. 8, 
1.0, 1.2, 1.6, 2.0, 2.4, and 2.7. 


Sea-Level Takeoff 

For sea -level takeoff with maximum afterburning, the nozzle total temperature was 
1945 K (3500° R), the total pressure was 3. 03 atmospheres (44. 5 psia), the primary flow 
rate was 287 kilograms per second (633 lbm/sec), and the nozzle pressure ratio Pq/Pq 
was 3. 03. 
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Hot -gas pressure distribution. - To calculate local velocity and subsequently the 
local heat flux, the plug surface pressure distribution is required. The pressure distri- 
bution on a 10° half -angle plug, operating at a nozzle pressure ratio of about 3. 0, is char- 
acterized by a series of over -expansions and recompressions along the plug surface. 
Figure 7 shows experimental pressures from three separate experimental studies. The 
Chenoweth data (ref. 6) were obtained on a model with a 12. 24-centimeter (4. 82-in. ) plug 
diameter and a full 100 -percent plug in a static test facility. The maximum gas total 
temperature in the study was 555 K (1460° R), and, for the data shown, a film-cooling slot 
was located at the 10-percent point on the plug. Reference 6 also showed that the position 
of the coolant slot did not significantly affect the plug pressure distribution. There was 
no secondary air used during the study. The Huntley data (ref. 4) were obtained on an un- 
cooled plug installed on a turbojet engine. The plug model was 43. 9 centimeters (17. 3 
in. ) in diameter, and, for the data shown, had a corrected secondary weight flow ratio of 
0. 071. The engine was run at rated speed with no afterburning, which resulted in 
primary-gas total temperature of 922 K (1660° R). The Bresnahan data (ref. 7) were 
also obtained on a 12. 24 -centimeter (4. 82 -in.) plug diameter model in both a static test 
facility and in a wind tunnel. Room -temperature air was used for both primary and sec- 
ondary streams; there was no secondary flow for the data shown. 


Nozzle Ref. 
pressure 
ratio, 

P 8^P0 

O 3.1 7 



Figure 7. - Comparison of predicted and experimental pressures on plug surface. Sea-level 
takeoff, Mach 0. 
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The data shown in figure 7 indicate that the plug static pressure distribution is a 
strong function of nozzle pressure ratio; this is also shown in the other data of refer- 
ences 4, 6, and 7. No general method has been found in the literature to calculate these 
pressure distributions. The pressure distribution is a function of the total temperatures 
of the primary stream and the secondary stream and a function of the total pressures of 
the two streams. The geometry, especially at the sonic point, is also critical. Finally, 
a mixing model must be assumed between the primary and secondary stream and between 
the secondary stream and the atmosphere in order to calculate the pressure profiles. No 
attempt was made in this study to solve this complex problem. The solid line in figure 7 
indicates the assumed pressure distribution used for this study. It was determined by the 
area-ratio method described in appendix B. A study was made of the effect of this as- 
sumption on the resulting wall temperatures and is discussed in the next section. 

Comparison of heat -transfer coefficients . - Heat -transfer coefficients must be known 
to calculate wall temperatures. Figure 8 shows a comparison of several methods of cal- 
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Figure 8. - Comparison of convective heat -transfer coefficients for sea-level takeoff and maximum afterburning. 


culating heat -transfer coefficients on the sting support and plug surface. The solid line 
was found using the fully developed pipe flow equation (ref. 10) and the area-ratio method 
of calculating pressure distribution. The short -dashed line represents the results of the 
Bartz boundary -layer analysis (ref. 11), again, using the area -ratio method to calculate 
pressures. The Bartz method was also run using measured pressures on the supersonic 
side, and the results are plotted as the dashed line in figure 8. The Bresnahan data 
shown in figure 7 were used for the pressure distribution. The effect of the overexpan- 
sions and recompressions is clearly seen on the heat -transfer coefficient with variations 
as high as 30 percent from the constant pressure assumption. 
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The effect of using these different heat -transfer coefficients on wall temperature is 
seen in figure 9. The sharp temperature gradients on the subsonic side of the plug result 
from the changes in the number of fins along the passage. These spikes would be 
smoothed in the actual case by conduction in the wall, which was neglected in this study. 
Neglecting wall conduction should not change the average wall temperature. Comparison 
of the two curves using the Bartz coefficients indicates that the difference between the 
constant -pressure assumption and the variable -pressure assumption results in a wall- 
temperature difference of only about 56 K (100° R) - or about 4 percent. The pipe flow 
equation yields the highest coefficients that occur at the nozzle throat, and, therefore, 
yields the highest, most conservative, wall tem per atures. For this reason, and since 
the pipe flow equation was much easier to work with, it was used to perform most of the 
calculations described in this report. 

The results of the turbulent flat -plate equation shown in figure 8 were used to calcu- 
late the convective heat transfer to the film -cooled extension. The results from the flat- 
plate equation are in reasonably good agreement with the coefficients used on the convec- 
tively cooled portion of the plug, as expected. 

Hot -gas temperature. - Another factor affecting the wall temperature is the radial 
hot -gas total -temperature profile. Figure 24 of reference 3 presents several gas- 
temperature profiles that were measured just downstream of the primary nozzle throat in 
the J-85-GE-13 engine plug nozzle tests. The profiles were similar in that each exhib- 
ited hotter gas temperatures at the center of the stream than near the plug surface or the 
primary nozzle lip. This was not a boundary -layer phenomenon, however, but was at- 
tributed to the annular burning pattern in the afterburner. As a result, some uncertainty 
exists as to which temperature to use as the driving temperature in the calculation of the 


Pipe flow equation; area-ratio 

pressure assumption 

Bartz; araa-satio pressure 

assumption 

Bartz; varying supersonic 

pressures 



Figure 9. - Effect of heat-transfer coefficient on wall temperature. Sea -level takeoff; maximum afterburn- 
ing; primary static temperature, 1945 K (3500^ R); 2 percent coolant from compressor discharge. 
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heat transfer to the plug wall. It was decided to use both the average or bulk gas temper- 
ature and an estimate of the gas temperature near the plug wall, for this study. The 
measured gas temperatures near the plug wall, shown in figure 24 of reference 3, were 
cross plotted as a function of average gas temperature and a linear relation resulted. 
Extrapolation to an average total temperature of 1945 K (3500° R) resulted in a gas total 
temperature near the wall of about 1667 K (3000° R). Each of these two gas temperatures 
were used as the driving temperature in the calculations for this report. 

Effect of coolant flow rate on wall temperatures . - Figure 10 shows the effect of 
varying the amount of cooling air on the plug wall temperatures. For this figure the pipe- 
flow equation was used to calculate the heat -transfer coefficients, and the area-ratio 
method was used to calculate local pressures and velocities. The air coolant for this fig- 
ure was assumed to have been obtained after the last stage of the compressor (tempera- 
ture, 667 K (1200° R), maximum pressure, 10. 2 atm (150 psia)). In figure 10(a), it can 
be seen that less than 2 -percent coolant is required to cool the plug below the maximum 
design wall temperature of 1220 K (2200° R), for a driving gas temperature of 1667 K 
(3000° R). In figure 10(b), for a driving gas of 1945 K (3500° R), about 3 -percent coolant 
would be required to maintain wall temperatures below 1220 K (2200° R). Note that com- 
pressor discharge air was assumed for the coolant; if the pressure requirements are sig- 
nificantly lower than the available 10. 2 atmospheres (150 psia), the air may be removed 
from an earlier stage on the compressor. The air temperature would then be reduced 
and the resulting wall temperatures could be reduced. 

Coolant pressure requirements . - Figure 11 presents the calculated pressure re- 
quirements as a function of coolant flow rate and location of film cooling slot (i. e. , the 
end of the convective-cooling passage and the beginning of the film-cooled extension). 
Again, the pipe flow equation was used to calculate the hot -gas heat transfer to the plug, 
and the area -ratio method was used to calculate pressure. The coolant was assumed to 
have been obtained from the compressor discharge. Figure 11(a) shows the calculated 
coolant total pressure at the end of the convective -cooling passages based on flow rate, 
area, and temperature. No data are shown for the 10 -percent slot location for flow rates 
greater than 2. 5 percent of the primary flow because the cooling air choked in the cooling 
passage, not at the end of the passage. 

Figure 11(b) shows the calculated total pressure required at the inlet to the coolant 
channels. Minor losses between fin sections were shown to be small in the data of refer- 
ence 3 and were, therefore, neglected. Also shown in the figure are the compressor - 
stage exit pressures. The results show that, if the plug were to be convectively cooled 
to the 75 percent point on the plug, cooling air would have to be obtained after the eighth 
stage of the compressor (3. 5 percent coolant). It was decided, however, to incorporate 
some film-cooled portion of the plug in order to (1) recover some thrust from the heated 
coolant, (2) reduce the plug weight, and (3) reduce the plug pressure drop. For a coolant 
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(a) Exit pressure. <b> Inlet pressure. 


Figure 11. - Effect of coolant flow rate and slot location on coolant total pressure. Cooling air assumed to be 
taken from the end of the compressor. Sea-level takeoff; maximum afterburning. 


discharge station at the 50 -percent point on the plug, a compressor discharge port after 
the fourth stage would provide cooling airflow up to about 3. 5 percent of the primary air- 
flow. Similarly, a discharge port after the third stage would provide up to about 3 per- 
cent of the primary airflow; figure 10 indicated that 3 -percent coolant would be the max- 
imum required for the convectively cooled portion of the plug (Tp = 1945 K (3500° R)) and 
that the actual coolant required would probably be less than 2 percent (T p = 1667 K 
(3000° R)). 

If the end of the convective-cooling passages could be located at the 25-percent point 
on the plug, a discharge port after the second stage of the compressor would supply about 
3-percent coolant to the plug. Similarly, for a plug with convection cooling to the 10 per- 
cent point and film cooling to the 75 -percent truncation point, a compressor discharge 
port after the first stage of the compressor would supply about 2. 5 percent of the engine 
primary airflow to the plug. 

For coolant flow rates greater than 2. 5 percent and the 10 -percent slot location, the 
coolant chokes upstream of the convective-cooling -passage exit. In other words, the 
minimum area is no longer the exit, but the junction between the sting support and plug. 
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coolant 



Figure 12. - Total temperature of coolant discharging from plug 
cooling channels. Sea -level takeoff; maximum afterburning 
cooling air obtained from end of compressor. 


For this reason these portions of the curves are not shown in figure 11. 

Temperature of coolant at film -cooling slot . - Let us now briefly discuss some of the 
cooling aspects of the film-cooled extension. Figure 12 presents calculated coolant total 
temperatures at the end of the convective -cooling passage. As in figure 10, the coolant 
was assumed to have been obtained from the end of the compressor (T = 667 K (1200° R)). 
These temperatures are required input for the film-cooling calculations and would be 
lower if cooling air can be taken from an earlier stage on the compressor, where the tem- 
perature is much lower. 

Effect of coolant temperature and flow rate on film cooling . - Figure 13 shows the ef- 
fect of cooling air temperature and flow rate on plug extension wall temperatures. The 
curves presented were calculated using hot -gas temperatures of 1945 K (3500° R) and 
1667 K (3000° R), a nozzle pressure ratio of 3.0, and a film-cooled plug extension from 
the 10 -percent point on the plug to the 75 -percent truncation point. Only those portions 
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Distance downstream of nozzle throat, percent 

(a) Primary total temperature, 1945 K (3500° R). (b) Primary total temperature, 1667 K (3000° R). 

Figure 13. - Effect of coolant total temperature and coolant flow rate on plug extension wall temperature. 10-Percent slot; wall emissivity, 0.65. Curves 
are shown within range of data correlation (ref. 6). 


of the curves within the range of the film-effectiveness data correlation of reference 6 
are shown. 

For a 6-percent coolant flow rate (fig. 13(a)) the plug wall temperature is approx- 
imately the same at the end of the plug, whether the coolant total temperature was 417 K 
(750° R) or 1111 K (2000° R). For the 2. 5- and 1-percent coolant flow-rate cases, the 
maximum design wall temperature will apparently be exceeded. Thus, it appears that 
lowering the coolant total temperature will have little effect on plug extension maximum 
wall temperatures for this application, except in the immediate vicinity of the film- 
injection slot. The effect of an increase in coolant flow rate, however, is characterized 
by a shifting of the wall -temperature curves downward and to the right. It appears from 
figure 13(a) that about 6 percent of the engine primary flow is required to film cool an ex- 
tension from the 10-percent point to the 75-percent truncation point using cooling air ob- 
tained at the end of the compressor and assuming that the total gas temperature is 1945 K 
(3500° R). Six percent cooling air is too high to maintain high cycle efficiency; therefore, 
using the 10-percent slot is impractical for these conditions. Similarly, a slot further 
upstream would require even more cooling air and for that reason was not considered in 
this study. 

If the effective primary total temperature is assumed to be 1667 K (3000° R), then 
the use of the 10-percent coolant slot position does appear possible. Figure 13(b) shows 
that the wall temperature will be approximately equal to the design wall temperature. 

Most of the wall cooling is being accomplished by wall -to -ambient radiation. The symbol 
shown in figure 13 indicates the wall temperature that would result with no film-cooling 
air (w /w = 0) ; all cooling would be accomplished by radiation. 

c p 
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Effect of slot location on wall temperatures . - Figure 14 shows the effect on wall 
temperature on film-cooling slot location using cooling air obtained at the end of the com- 
pressor. Figure 14(a) shows the results of convection cooling to the 50-percent point on 
the plug and film cooling to the 75 -percent truncation point. Figure 14(b) shows the re- 
sults of convection cooling to the 25 -percent point, and figure 14(c) shows the results for 
convection cooling to the 10-percent point. Two and one-half percent of the primary air 
was assumed for this figure, and the pipe flow equation was used to calculate the heat- 
transfer coefficients to the convectively cooled portion of the plug. Figure 14 shows, for 
the convectively cooled portion, that an effective gas temperature of 1945 K (3500° R) 
results in a marginal design with a maximum wall temperature of about 1258 K (2265° R) 
at the nozzle throat. Using an effective gas temperature of 1667 K (3000° R), which is 
probably more realistic, reduces the nozzle throat wall temperature to about 1111 K 
(2000° R). 

The film-cooled extension temperatures for the case with film cooling from the 50- 
percent point are well below the design temperature limit of 1222 K (2200° R). The effect 
of wall emissivity e w is small over the expected range for the materials involved. The 
film -cooled wall temperatures for the slot located at the 25-percent point are shown in 
figure 14(b). These temperatures are marginal for the 1945 K (3500° R) gas temperature 
but adequate for the 1667 K (3000° R) gas temperature. Similarly, the wall temperatures 
with convection cooling to the 10-percent point on the plug (fig. 14(c)) appear to be too 
high for the higher gas temperature but reasonable for the lower gas temperature. 

Now, consider what would happen to these extension temperatures if lower tempera- 
ture cooling air could be used to cool the plug (from an earlier compressor stage, for 
example). It was shown in figure 13 that cooiant-inlet temperature had little effect on 
maximum extension temperature for a coolant slot located at the 10-percent point because 
the downstream film effectiveness was small and most of the cooling was accomplished by 
radiation. Therefore, reducing the coolant -inlet temperature would not significantly re- 
duce the maximum -extension wall temperature for the 10-percent slot case. The same 
reasoning applies for the slot at the 25-percent point. On the other hand, the configura- 
tion with the extension located at the 50-percent point would be favorably affected by a re- 
duction in coolant -inlet temperature because of the short length of the film -cooled portion. 

Resulting wall temperatures using interstage compressor bleed . - The cooling re- 
quirements for the sea -level -takeoff flight condition are summarized in figure 15. The 
experimental pressure distribution of Bresnahan (ref. 7) was used on the supersonic side 
of the plug to calculate local Mach numbers. Then the Bartz boundary -layer method was 
used to calculate gas -to -wall convection heat transfer. These assumptions are believed 
to be the most realistic for final calculations. The curves shown were calculated using 
the higher gas temperature, 1945 K (3500° R). Obviously, lower wall temperatures 
would result if the lower gas temperature were used. The lower curve was calculated 
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Figure 15. - Resulting wall temperatures using compressor Infer stage bleed. Sea -level takeoff; maximum 
afterburning, primary total temperature, 1945 K (3500^ R). 


using 2. 5 percent of the primary airflow to cool the plug and sting support, removed after 
the third stage of a nine-stage compressor. The temperature of the cooling air at this 
stage is only 406 K (731° R), and the pressure is 4. 60 atmospheres (67. 6 psia). The cal- 
culated total pressure required at the inlet to the sting support is only about 3. 28 atmos- 
pheres (48. 2 psia); therefore, 1. 32 atmospheres (19. 4 psia) are available to duct the air 
from the compressor to the sting support inlet. Since these wall temperatures are well 
below the maximum wall temperature, this configuration is feasible. 

The upper curve in figure 15 was calculated using 2 percent of the primary airflow 
removed after the second stage of the compressor. The temperature and pressure of the 
air available at this point are 367 K (660° R) and 3.4 atmospheres (50 psia), respectively. 
The calculated pressure required at the inlet to the sting support for this configuration is 
only 2. 72 atmospheres (40.0 psia). Again, 0. 68 atmosphere (10. psia) is available be- 
tween the compressor bleed port and the sting support inlet. These wall temperatures 
are also well below the maximum wall temperature, so this configuration also appears to 
be practical. The film-cooled extension temperatures are also satisfactory. Since there 
are performance gains to be made from removing this air from as early a stage as possi- 
ble, the latter configuration would be favored. 


Supersonic Cruise 

Upon selecting a cooling configuration for one flight condition (takeoff, in this case), 
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it becomes necessary to evaluate the cooling requirements at other flight conditions. The 
supersonic -cruise condition assumed for this study included a flight Mach number of 2.7, 
a primary flow rate of 186 kilograms per second (410 lbm/sec), a nozzle total pressure 
of 1. 97 atmospheres (29 psia) at beginning of cruise and 1. 35 atmospheres (20 psia) at 
end of cruise, and a nozzle pressure ratio of 27.6. The afterburner was presumed to be 
on at a low power setting, resulting in a nozzle total temperature that varied from 1347 
to 1286 K (2425° to 2314° R) from the start of cruise until the end of cruise, respectively. 

o 

The primary nozzle flow area was reduced from 10 600 square centimeters (1642 in. ) 

o 

for maximum afterburning to 8500 square centimeters (1317 in. ) at supersonic cruise. 

Plug static pressure distribution . - Figure 16 shows measured plug static pressures 
for 10° half -angle conical plugs operating at high pressure ratios. The Bresnahan data 
(ref. 7) were obtained on a 12. 24 -centimeter (4. 82 -in. ) plug diameter model in both a 
wind tunnel and in a static -test facility. The location shown for the end of the cylindrical 
ejector for these data is approximately optimum for this nozzle pressure ratio. No sec- 
ondary flow was used for the data shown, however. The Huntley data (ref. 4) were ob- 
tained from an uncooled plug attached to an afterburning turbojet; the engine could only be 
operated in the nonafterburning mode, however, because of the temperature limit on the 
uncooled plug. The corrected secondary weight flow ratio was 6. 3 percent for these data, 
and the secondary shroud location was not optimized. The Clark data (ref. 3) were ob- 
tained on a cooled plug with nonafterburning for the case shown. The cylindrical ejector 
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Figure 16. - Comparison of predicted and experimental pressures on plug surface. Super- 
sonic cruise; Mach 2. 7. 
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position was not optimized for the nozzle pressure ratio. Secondary air was used, 
however. 

The calculated pressure distribution, using the approximate, area-ratio method, 
yields reasonable results up to about the 10-percent point on the plug and between the 30- 
percent point and the 50-percent point. From about the 10-percent point to the 30 percent 
point, the estimate is too high, by a factor of about 2 at the 25-percent point. 

Comparison of predicted heat-transfer coefficients. - Figure 17 compares the var- 
ious methods used to predict heat -transfer coefficients for the supersonic -cruise flight 
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Figure 17. - Comparison of heat-transfer coefficients. Supersonic cruise; primary total temperature, 1345 K 
(2425° R); Mach 2. 7 : coolant total temperature, 533 K (960° R). 


regime. The dashed line in figure 17 shows the coefficients using the pipe flow equation 
and pressures using the area-ratio method. The solid line on figure 17 was found by us- 
ing the Bartz method, again using the area-ratio method of predicting pressures. The 
short -dashed line in the figure was calculated using the Bartz method and the experi- 
mental pressures of Bresnahan (see fig. 16) to evaluate the error introduced by using the 
area-ratio method of calculating pressures. Comparison of the two Bartz curves shows 
that the heat -transfer coefficient distribution is similar to the pressure distribution as- 
sumed. Also, the pipe flow equation predicts a higher heat -transfer coefficient over the 
entire plug. Therefore, the use of the pipe-flow equation should yield a conservative 
design. 

Effect of hot -gas temperature profile. - A hot -gas temperature profile will exist in 
the primary gas stream even for the low level of afterburning assumed (see ref. 3). The 
relation found from reference 3 between average gas temperature and gas temperature 
near the wall yields, for an average gas temperature of 1347 K (2425° R), a gas temper- 
ature near the plug of about 1166 K (2100° R). Thus, if the afterburner burner profile 
found on the J-85-GE-13 afterburning turbojet of reference 3 can be duplicated on the 
larger supersonic -cruise aircraft, no plug cooling will be required for supersonic cruise, 
and plug wall temperatures will be about 1166 K (2100° R). In the discussion that follows 
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an average gas temperature of 1347 K (2425° R) is assumed to be the driving temperature, 
and the cooling air is assumed to have been obtained from the engine inlet with a ram 
pressure of 1. 52 atmospheres (22. 3 psia) and temperature of 533 K (960° R). 

Ram-air coolant requirements. - Figure 18 summarizes the effect of ram-air cool- 
ant flow rate on both plug wall temperature at the primary throat and calculated plug- 
inlet coolant pressure. The abscissa shows the coolant flow rate as a percent of the en- 
gine primary airflow and the right ordinate is a measure of the required inlet pressure 
at the cooling passage inlet station. The left ordinate shows plug wall temperature at the 
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Figure 18. - Effect of coolant flow rate on plug throat wall temperature and required coolant 
inlet pressure. Supersonic cruise point; primary total temperature, 1345 K (2425° R). 
Mach 2. 7; ram cooling air ; coolant total temperature, 533 K (960° R). 


nozzle throat (station 8). The calculations were made using the area-ratio method to cal- 
culate pressures on the supersonic side; heat-transfer coefficients were calculated using 
the pipe flow equation. 

The figure indicates that about 1. 5 percent of the primary airflow may be obtained 
using a ram inlet and that this flow would result in a plug wall temperature at the primary 
throat of about 978 K (1760° R). Restricting the flow would result in higher wall temper- 
atures, but a flow rate as low as 0. 5 percent results in a temperature of only about 
1128 K (2030° R) - well below the design wall temperature of 1222 K (2200° R). Obvi- 
ously, the selection of the correct flow rate for this flight condition will also include a 
life criteria; lower temperatures generally result in longer lives. Operation with the 
coolant flow rate between 0. 5 and 1. 5 percent of primary flow would result in plug tem- 
peratures between the two curves in figure 19, which should yield satisfactory life. 
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Figure 19. - Wall temperatures at supersonic cruise. Primary total temperature, 1345 K (2425° R); Mach 2. 7; ram air 
coolant; pipe flow equation. 


Maximum Acceleration 

Figure 20 compares plug wall temperature at the nozzle primary throat for the entire 
Mach number range up to 2.7 for maximum afterburning. Two -percent coolant was as- 
sumed to have been taken from the engine after the second stage of the compressor. For 
these calculations, local velocity and pressure were found using the area-ratio method, 
and heat -transfer coefficients were calculated using the pipe flow equation. The figure 
shows that the maximum plug wall temperature T w g will be below the design temper- 
ature except for takeoff and above Mach 1. 85 for the 1945 K (3500° R) curve. When the 
gas temperature profile is taken into account (the 1667 K (3000° R) curve in fig. 20), the 
maximum plug wall temperature is predicted to be below 1162 K (2093° R). Actually, on 
a typical mission, the pilot would throttle the afterburner back as he approached Mach 
2. 7, thus reducing the plug wall temperature. It appears that 2-percent coolant, re- 
moved from the compressor after the second stage, should provide adequate plug cooling 
during maximum afterburning. 
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Figure 20. - Plug throat wall temperature as function of Mach number. Maximum 
afterburning; 2 percent coolant from compressor second stage. 


OTHER CONSIDERATIONS 

As discussed in the INTRODUCTION of this report, certain other considerations are 
important in a nozzle design. The plug weight, assuming convection cooling to the 50 per- 
cent point, was calculated by considering the outer wall, fins, and inner wall and was 
found to be about 616 kilograms (1361 lbm). For four engines this is about 0. 7 percent of 
the entire aircraft weight. Breaking this weight down into the various components: 


Plug and sting support outer wall, kg (lbm) 219 (482) 

Plug and sting support inner wall, kg (lbm) 56 (123) 

Plug and sting support fins, kg (lbm) 315 (694) 

Upstream of primary throat, kg (lbm) 81 (178) 

Downstream of primary throat, kg (lbm) 219 (482) 

Sting support, kg (lbm) 15 (34) 

Film -cooled extension, kg (lbm) 26 (57) 

Total, kg (lbm) 616 (1361) 


Notice that about 35 percent of the weight is in the fins downstream of the primary 
throat. Thus, if the convective cooling passage ended at the 25 -percent point and the plug 
extension would film cool to the 75-percent truncation point, about 100 kilograms 
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(221 lbm) could be saved per engine. Similarly, if the extension would film cool all the 
way from the 10-percent point to the 75-percent point, about 186 kilograms (410 lbm) 
could be saved per engine. 

Lateral vibration was checked using the method of reference 12. The method as- 
sumes that the first natural frequency in the lateral mode can be found by assuming that 
the sting support extends to the center of gravity of the plug and that the mass of the plug 
is concentrated at the center of gravity. Also, 23. 6 percent of the weight of the sting 
support is assumed to be concentrated at the plug center of gravity as outlined in the 

method of reference 12. Young's modulus was assumed to be 1. 59X10 11 newtons per 

£» 

square meter (23. 1x10 psi), and the moment of inertia of the sting, including the fins, 
was calculated to be 21 500 cm^ (517 in.^). Thus, the first natural frequency was calcu- 
lated to be 6. 65 hertz. 


Stresses in the support tube wall were checked using elementary strength of mate- 

5 

rials theories. The longitudinal stress caused by an internal pressure of 3.45x10 new- 



7 

circumferential stress resulting from the same internal pressure is 2. 76x10 newtons 

per square meter (4000 psi). Bending stress was calculated by assuming an equivalent 

static load at the plug center of gravity resulting from an acceleration of four gravities. 

The maximum bending stress at the junction with the turbine frame was found to be 
7 

5. 86x10 newtons per square meter (8500 psi). This stress combines with the longitu- 

7 

dinal stress to yield a maximum tensile stress of 7. 24x10 newtons per square meter 
(10 500 psi). Thus, to maintain a 1000 -hour rupture strength, the temperature of the 
wall at this point must be kept below about 1110 K (2000° R). Similarly, the bending 
stress was checked at the junction between the sting support and the plug. The bending 
stress was found to be 3.38x10 newtons per square meter (4900 psi), and the total ten- 
sile stress 4. 76x10^ newtons per square meter (6900 psi). Based on the calculated wall 
temperature of about 1166 K (2100° R) (see fig. 15) at this point, the rupture life should 
be between 100 and 1000 hours. This life can be improved by improving the cooling in 
this area or by building stiffeners into the sting support. 

Buckling of the support tube wall was also checked and was not found to be a problem. 
The purpose of these preliminary stress calculations was to determine whether or 
not the wall size selected was reasonable. The numbers indicate that it is. 

SUMMARY OF RESULTS 

The important results of this thermal design study may be summarized as follows : 

1. An air-cooled plug-nozzle system for an afterburning turbojet engine for a 
supersonic -cruise aircraft appears to be practical using low-cycle-penalty cooling air. 
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During maximum afterburning (i. e. , sea-level takeoff and acceleration up to Mach 2. 7), 

2 percent of the primary airflow may be removed after the second stage of a nine-stage 
compressor and should result in plug wall temperatures below 1222 K (2200° R). 

2. During supersonic -cruise operation, ram air may be used to cool the plug and 
sting support. About 1. 5 percent of the primary airflow may be obtained at ram pres- 
sure and temperatures , which would result in maximum plug wall temperatures of about 
978 K (1760° R), when it is assumed that the average gas temperature is the driving tem- 
perature. Restricting the coolant flow to about 0. 5 percent would result in maximum 
wall temperatures of about 1133 K (2040° R). If the hot-gas temperature profile can be 
tailored so that the gas temperature near the plug wall is below 1222 K (2200° R), then no 
cooling air should be required during supersonic cruise. 

3. The configuration selected includes a 10° half -angle plug that is sting supported to 
the turbine frame. Convection cooling on the sting support and plug up to the 50 -percent 
point on the plug (half-way between the primary nozzle and the end of a full plug) should 
provide adequate cooling with low pressure drop. At the 50 -percent point, the air cool- 
ant would be injected as a film over a 10° half -angle extension that is truncated at the 75- 
percent point. Thus, some thrust may be recovered from the heated coolant. 

4. Preliminary stress calculations indicate that the wall size selected is reasonable. 
The plug weight was calculated to be about 616 kilograms (1361 lbm) for a plug with con- 
vection cooling to the 50 -percent point and film cooling to the 75 -per cent point. Lateral 
vibration of the plug was checked for this configuration, and the first natural frequency of 
the plug was found to be 6. 65 hertz. 

5. A configuration that includes convection cooling to the 10-percent point on the plug 
and film cooling to the 75-percent point also appears practical if the afterburner temper- 
ature profile can be tailored so that the driving temperature near the plug wall can be 
maintained at 1667 K (3000° R). This configuration would have a 186 kilograms (410 lbm) 
weight advantage over the configuration with convection cooling to the 50-percent point. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, October 14, 1971, 

764-74. 
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APPENDIX A 


SYMBOLS 


A 

area 

P 

total pressure 

* 

A 

critical area 

Pr 

Prandtl number 

a 

absorptivity 

P 

static pressure 

a gw 

absorptivity of gas from a source 

T 

P 

partial pressure 


at T 

w 

q 

heat-transfer rate per unit area 

C 

correction factor in radiation 

R 

gas constant 


(eq. (C6)) 




Re 

Reynolds number 

C D 

discharge coefficient, 



Aetuabflow/Ideal flow 

R i_ 

inner radius of hot -gas flow 



channel 

c p 

specific heat at constant pressure 

R n 

outer radius of hot -gas flow 

D 

diameter 

0 

channel 

D h 

hydraulic diameter 

s 

film-cooling step height 

F 

radiation configuration factor 

Sp 

fin spacing 

f 

friction factor 

T 

total temperature 

f/a 

fuel -air ratio 

T 

w 

wall temperature 

G 

specific flow rate, w/A 

t’ 

w 

adiabatic wall temperature 

g 

Newton's conversion constant 

T 

rec 

recovery temperature 

H 

fin height 

t 

static temperature 

H* 

length of equipotential line 

V 

velocity 


(i. e. , passage height) 




w 

fin base width 

h 

heat -transfer coefficient 

w 

mass flow rate 

I 

station index 

X 

distance along plug surface 

J 

Joule's constant 


downstream of slot 

k 

thermal conductivity 

X 

distance from base of fin 

L 

length 

a 

thermal diffusivity 

h> 

mean beam path length 

y 

ratio of specific heats 

M 

Mach number 

£ 

emissivity 
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A 

recovery factor 

i 

inside 

M 

dynamic viscosity 

in 

inlet to convective cooling slot 

P 

density 

m 

momentum 

a 

Stefan-Boltzman constant 

max 

maximum 

T 

thickness 

out 

outside 

Subscripts: 

P 

primary 

am 

ambient radiation conditions 

Pi 

plug 

av 

average 

r 

radiation 

c 

coolant 

ref 

reference properties 

CO 

carbon dioxide 

si 

coolant slot 

conv 

convection 

w 

wall 

eff 

effective 

wa 

water vapor 

exit 

convective passage exit 
(coolant slot) 

X 

distance along plug surface down- 
stream of slot 

f 

film 

1 

hot -gas -to -wall radiation 

fl 

flame 

2 

wall -to -coolant forced convection 

fin 

fin 

3 

wall -to -ambient radiation 

fr 

friction 

8 

primary nozzle station 

I 

station index 

0 

ambient flow conditions 
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appendix b 


CONVECTION HEAT -TRANSFER AND PRESSURE-DROP ANALYSIS 

The procedure used to analyze the plug -nozzle system was developed into two com- 
putor programs (FORTRAN IV): One program was used for the convection -cooled por- 
tions of the sting support and plug, and the second program analyzed the film-cooled plug 
extension (appendix C). The program used for the convection -cooling analysis assumes 
one -dimensional, steady-state heat flow through the outer wall and fins into the coolant. 
Transport properties of the hot gas were obtained from reference 13. A fuel-air ratio of 
0. 0566 was used for the maximum afterburning cases and 0. 02 was used for supersonic 
cruise. 


Input Data 

The heat-transfer model used for the convection -cooled portion of the sting support 
and plug is shown in figure 6(a). The thermal conductivity of the outer wall and fins were 
assumed to be constant: 

k w = 32.9 J/(sec)(m)(K) (19. Btu/(hr)(ft)(°R)) 

for a typical high strength nickel -base alloy and 

k fin = 55.3 J/(sec)(m)(K) (32. Btu/(hr)(ft)(°R)) 

for nickel. The cooling passages were divided into equal 5. 08-centimeter (2. 0-in. ) in- 
crements, and station calculations were performed on each increment. Coolant inlet 
total temperature and flow rate, hot -gas total conditions, and plug geometry were input 
to the program. 


Heat Flux Loop 

An energy balance was performed at each station that resulted in an estimate of the 
local wall temperature. 

Coolant side . - At each station the following parameters were calculated: (1) Fin 
spacing, (2) coolant flow area, (3) coolant flow perimeter, (4) hydraulic diameter, 

(5) Reynolds number, and (6) heat -transfer coefficient. The coolant transport properties 
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were approximated from air tables as a function of the coolant total temperature. The 
coolant flow was turbulent for all conditions studied. The Dittus -Boelter correlation for 
turbulent flow (ref. 14) was used to calculate heat -transfer coefficients 

h c,i = 5 ^-(°' 023 Re 2 ;i 8pr °' 4 ) (B1) 

h,I 

The effect of the high -conductivity fins was included by solving the one -dimensional con- 
duction equation with the following boundary conditions (ref. 15): 

(1) The temperature at the base of the fin is constant and equal to T .. 

(2) The temperature gradient at the end of the fins is zero (dT/dx = 0). 

(3) The heat -transfer coefficient is constant and the same on the inner wall and fins. 
The effective heat -transfer coefficient on the coolant side can be shown to be 



Hot -gas side. - On the hot-gas side the local calculated area ratios were used to 
calculate local hot -gas conditions assuming a specific heat ratio of 1.3 (ref. 16). In the 
subsonic region of the nozzle, this method was believed to yield a good prediction of hot- 
gas static pressure, static temperature, local velocity, and density. Problems develop 
in the supersonic region, however. The flow and, hence, pressure field in the super- 
sonic region of the nozzle is a function of 

(1) Nozzle pressure ratio 

(2) Secondary shroud position 

(3) Secondary airflow rate 

(4) Primary gas flow rate 

(5) Temperature of primary and secondary streams. 

A method has been developed at Lewis to evaluate the performance of a fully ex- 
panded plug nozzle with secondary flow (ref. 1). A comparison of static pressures - 
experimental and predicted - is shown in reference 3. The program was developed for 
very high pressure ratios, that is, fully expanded flow. There is no known method for 
accurately predicting pressure distributions (and, hence, local velocities) for the lower 
pressure ratio cases. Therefore, a method was required for this study to evaluate the 
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lower pressure ratio cases (nozzle pressure ratio is shown in fig. 1(b) as a function of 
flight Mach number and varies from 3 to about 27). 

The prediction method selected for this study yields, at best, a rather crude approx- 
imation of the local pressure distribution. The technique requires a guess at the flow 
boundary between the primary and secondary streams. Specifically, a drawing was made 
of the plug, primary nozzle lip, and secondary shroud. Then, the flow boundary between 
the primary stream and the secondary stream was estimated and sketched on the figure. 
Figure 21 shows the sketch for this nozzle. Equipotential lines were then estimated and 
drawn from the plug stations to the outer flow boundary. The secondary shroud diameter 
was assumed to be the outer limit of the equipotential lines, even downstream of the end 
of the secondary shroud. Outer radius and passage height were then scaled from the 
drawing, and the local area ratio was calculated. 



* 

A 


tt(R. + R o )H’ 
A 8 


(B3) 


Figure 21(b) shows the equipotential lines estimated for the supersonic side of the nozzle. 
This estimate was simply a first approximation; a check was made of the calculated local 
static pressure to be sure it was above the back pressure. For the low -nozzle-pressure 
ratio cases, the calculated static pressure dropped below the back pressure at some sta- 
tion along the plug downstream of the primary nozzle lip. When this occurred, the local 
pressure was set equal to the back pressure, and the area ratio was set equal to the max- 
imum for this pressure ratio: 




(y+l)/2(y-l) 


(B4) 


where 


M = 1 / 2 

Wr-D/r _ - 

yr - l 



(B5) 


Then, the local outside radius, passage height, and hydraulic diameter were recalculated 
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Figure 21. - Equipotential lines used to calculate areas. 





with the assumption that the primary flow remains attached to the plug: 
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The effective gas temperature was then calculated to account for high-speed flow effects. 


T eff,I _t p,I + A(T 8 -t p,I* 


(B9) 


where 


A = VPTpgf for turbulent flow (BIO) 

The reference Prandtl number was calculated as a function of reference temperature 
(ref. 17). 


T ref, i ~ °- 5 T w , out, I + °- 28 *p, I + 0> 22 T eff , I 


(Bll) 


Finally, a simple equation, originally developed for fully developed pipe flow, was in- 
cluded to calculate convective heat -transfer coefficients from the primary gas to the plug 
wall (ref. 10). 


V I - 0. 026 (Re ref )°- 8 (Pr ref ) 0 ' 3 (B12) 

Although this equation was believed to be reasonably accurate and was very easy to use, 
it was decided to compare the results from this equation with a more sophisticated 
boundary -layer procedure. Therefore, an option was included in the program to read in 
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values of local heat -transfer coefficient. Then, any of several existing boundary -layer 
programs could be used to calculate gas-to-wall heat transfer. The method selected for 
this study was developed by Bartz (ref. 11). The Bartz method solves the integral mo- 
mentum and energy equations of the boundary layer in conjunction with a modified Von 
Karman momentum-heat analogy. The input to the Bartz program includes plug geometry 
and estimated wall temperature, total gas conditions, initial momentum boundary-layer 
thickness, and local Mach number. The problem that exists in using the Bartz program, 
as well as any other method of this type is the required local velocity or Mach number of 
the supersonic stream. As pointed out earlier, the estimate of local pressure (and, 
hence, velocity) is crude, at best, and may introduce more error than the difference be- 
tween the Bartz coefficients and the pipe-flow coefficients. 

The convective heat transfer from the gas to the wall is thus found: 

q conv = h p, I (T eff , I " T w, out, I ) (B13) 

Radiation. - Radiation from the hot combustion gases to the plug wall was calculated 
using a one -dimensional method described in reference 18. 


ct _ 1 / 1 . \ xl. 5/. 2. 5 rp 2 . 5 \ 

q r,fl ~ g CT 1 w^fl hi ‘ T w, out, 1/ 

The absorptivity of the wall a was assumed to be equal to 0. 80 for this study. 

w 

flame emissivity was found from 


(B14) 


The 


e fl = 1 - exp £-18. 5 P^ ^ T" 1 ' 5 j (B15) 

This empirical correlation was developed for a nonluminous flame, and the model 
was said to be accurate when axial temperature gradients and end effects are small and 
the length-to-diameter ratio of the flame tube is large. The mean beam path length was 
defined to be 




3, 6 (Radiating volume) 
(Receiving surface area) 


(B16) 


For an axial increment AL this suggests that the mean beam path length is approx- 
imately equal to the local hydraulic diameter. 




3, 6 (Radiating volume) 
(Receiving surface area) 


4 (Area) AL _ ^ 

— T 

(Perimeter) AL 


(B17) 
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Since the local hydraulic diameter varied considerably along the plug surface of the 
system studied and the mean beam path length is approximately equal to the hydraulic 
diameter, the local hydraulic diameter was taken to be the mean beam path length for 
this study. This is a good approximation along the sting support, where the flame radia- 
tion is highest. The flame radiation decreases considerably downstream, especially near 
the end of the supersonic plug surface. 

Radiation from the sting support and the convectively cooled portion of the plug to the 
afterburner liner, the primary nozzle, the secondary shroud, and the atmosphere were 
neglected for this study. The afterburner liner (fig. 2) and iris primary would be ap- 
proximately the same temperature as the sting support and plug surface; the net inter- 
change of radiant energy should therefore be negligible. The translating secondary 
shroud will be retracted for acceleration to about Mach 1. 2 and extended at higher Mach 
numbers. The temperature of the secondary shroud is a function of the amount of 
secondary air used, and this variable was not considered in this study. The net inter - 
change should be from t he plug to the secondary shroud, however . — This is al so the case 
for the radiation from the plug wall to the atmosphere. This term is only important near 
the end of the plug and, if included, would result in lower predicted plug wall temper- 
atures in this region. Since the throat region was the critical convection-cooling region 
for the cases studied, the inclusion of this radiation term would not have altered the con- 
clusions reached. 

Wall temperature. - The heat entering the plug wall locally, then, was taken to be 
the sum of the convection from the hot combustion gases and radiation from the flame: 


q_ q conv + q r,fl 

The inner wall temperature T z (see fig. 6) was then calculated: 

W, 1 


(B18) 


T = T 
w,i w,out 


r w fin 


-q — + 


i k w k fin/ 


(B19) 


This equation assumes that the fin base occupies the entire space between fins. This is a 
reasonable assumption since, for the geometry selected, the distance between fins 
(Sp - W) is only about 20 percent greater than the fin width W on the average. 

The coolant total temperature rise to the next station is, then, 


AT c 


_ qA 

WC P 


(B20) 
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Thus, the coolant total temperature at station I + 1 is calculated to be 


T c,I+l " T c,I + AT c 


(B21) 


The heat flux loop (eqs. (Bl) to (B21)) was repeated for each station from the coolant inlet 
to the end of the convective -cooling channel. 


Coolant Pressure-Drop Loop 

After completing the heat -flux calculations to obtain coolant temperatures, the cool- 
ant pressure calculations were performed from the discharge end of the passage, 
upstream, to the coolant inlet point. Choked flow is initially assumed at the end of the 
coolant channel. The total temperature, area, and flow rate are known. The static tem- 
perature was then found: 


t = 



And the total pressure was found from Fliegner's formula (ref. 19): 


(B22) 


P - 



0. 532 AC d 


(B23) 


Nonisentropic flow was accounted for by including a discharge coefficient C D of 0. 98. 
Next, the static pressure at the choked flow station was calculated: 


p = P |l + ^ ^ (B24) 

This static pressure was then compared with the nozzle back pressure Pq. If the calcu- 
lated static pressure was greater than p^, the flow was, in fact, choked. If not, the 
flow was not choked, and the passage exit conditions were found by setting the static 
pressure equal to the back pressure and calculating the local static temperature and total 
pressure. 
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(B25) 


P = p 
t 



(B26) 

(B27) 


The static temperature is initially set equal to the total temperature, and an iteration is 
required. Thus, the state of the fluid at the end of the convection -cooled channel is com- 
pletely defined. 

The following pressure-drop calculations were then made lor each station from the 
passage exit upstream to the coolant inlet station. Average Reynolds number, specific 
flow rate, hydraulic diameter, and station length were found and friction factor calculated 
assuming turbulent flow in smooth tubes (ref. 14). 


f = 0.0014 + 


0. 125 


Re 


0.32 

av 


(B28) 


Initially, the static pressure at station I is assumed to be equal to the static pressure 
at station I + 1. Then the static temperature at station I was calculated from the quad- 
ratic equation 


ms. 


I pfrgj 


tj + tj - Tj = o 


(B29) 


The static -pressure drop was found by summing the momentum and friction pressure 
drops: 


Ap 


m 


w 2 / 1 

A av g WI+1 A I+1 



(B30) 
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(B31) 


Ap {r = 



gD 


h,av 


4f AL 
( ^I + p I+l^ 


Thus, 


Pi = Pj +1 + Ap m + Ap fr (B32) 

An iteration was required until the calculated pj was equal to the assumed Pj. On con- 
vergence, the Mach number and total pressure were calculated explicitly: 


M 


I " 


Pj yrgRtj 


Expansion and contraction losses between fin sections were neglected as well as entrance 
losses at the turbine frame. 
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APPENDIX C 


FILM-COOLING HEAT -TRANSFER ANALYSIS 

The procedure used for predicting plug surface temperatures resulting from film 
cooling the aft end of the plug nozzle was based on the method of reference 8; it used the 
wall heat balance shown in figure 6(b). The internal plug temperature was assumed to be 
equal to the plug surface temperature. The axial conduction rate in the wall was as- 
sumed to be negligible. Also, the heat -transfer areas in the heat balance were assumed 
to be equal. 

A heat balance on the wall (fig. 6(b)) was assumed wherein the hot -gas -to -wall ra- 
diation q i was equal to the sum of the wall -to -coolant forced convection qg and the 
external -to -ambient radiation qg 


q t = q 2 + q 3 (Cl) 

The hot -gas -to -wall radiation q^ was calculated from the expression (ref. 14) 

q l-(^)(Vp,av-V T t) < C2 > 

The wall -to -coolant forced convection ^ was calculated from 

q 2 = h 2 (T w - t;> (CS) 

The external-to-ambient radiation q^ was calculated from 

is'vH-c) < c4 > 

By substituting the values for q^, qg, and q^ as defined in equations (C2), (C3), and 
(C4), respectively, into equation (Cl), an expression was obtained wherein all terms in 
the equation were determinable or expressed as a function of T . Equation (Cl) was 
then solved using the Newton-Raphson method of iteration to determine the wall temper- 
ature assinning an initial estimate of the adiabatic wall temperature T^ for T w (qg = 0 
for this assumption). The values used in equation (C2) to (C4) were determined as dis- 
cussed in the next section. 
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Internal Flow Properties 


The primary flow field was determined by assuming one-dimensional isentropic flow, 
hot-gas properties for the assumed ASTM A-l/air mixture, and uniform wall static pres- 
sure across the nozzle cross section. The wall static -pres sure profile was assumed in 
accordance with the calculated pressure distribution line (figs. 7 and 16). The axial dis- 
tance between the slot and the end of the plug was divided into equal increments. Local 
hot -gas flow properties at each station were determined and used in the heat balance 
(eqs. (Cl) to (C4)) as required. 

The coolant-flow field was determined by assuming that the air coming out of the 
plug convection -cooling passages flowed between the primary gas and the plug wall. To 
account for primary and coolant -stream thermal mixing, the local coolant properties 
were determined as follows. A semi -empirical film-cooling correlation was used to esti- 
mate a coolant -stream adiabatic wall temperature T^. This adiabatic wall temperature 
was assumed to represent the recovery temperature T rec c of the coolant stream. Us- 
ing the expression (ref. 20) 


T 

rec, c 



+ Pr 


1/3 

c 



2C gJ 
p c 


(C5) 


together with the equation of state, continuity, and the air property tables, the local 
coolant-stream temperature t £ at each station was determined. Local coolant proper- 
ties were determined at the local coolant temperature. The coolant -stream cross- 
sectional area was calculated by assuming isentropic flow from the coolant slot to match 
the wall static -pressure profile (figs. 7 and 16, calculated pressure distribution line). 

All coolant -slot exit properties were calculated using isentropic equations, continu- 
ity, air property tables, the coolant -slot static temperature t the total coolant tem- 
perature T , and the total coolant pressure P , which was obtained from the convection- 
cooling analysis discussed in appendix B. 


Hot-Gas-to-Wall Radiation 

The primary gases were assumed to result from the combustion of ASTM A-l (fuel) 
and air. Nonluminous radiation from the combustion products carbon dioxide and water 
vapor were considered and are accounted for in equation (C2). In order to solve equa- 
tion (C2) (taken from ref. 14), values for t p &v , e , a , and e w were required. The 

hot gas was assumed to radiate at each station at an average temperature t which 

p, av 

was obtained by averaging all local primary static temperatures from the slot to the end 
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of the plug. The values of emissivity e p and absorptivity a^ for hot-gas radiation 
were estimated as outlined in reference 14 and as shown by equations (C6) and (C7): 


I 

i 


V 


e co^co + Sva^wa 


Ae 


(C6) 


where 


e co = ^p’PcoV 
e wa " ^p’Pwa^b^ 

and where C CQ is the correction for total-pressure effect for carbon dioxide, C wa is 
the correction for total -pres sure effect for water vapor, and Ae is the correction for a 
mixture of COg and HgO. The absorptivity, 


where 


a = a 4- a 
gw co wa 


Aa 



(C7) 


a wa ^w’^wa^k 

and where Aa is the correction for a mixture COg and HgO. Using an average hot -gas 
static temperature and static pressure at the stations being considered, the fuel-air 
ratio, and the method of reference 21, the mole fractional values for the products of com- 
bustion (for assumed equilibrium composition) were determined. The average hot-gas 
static pressure was multiplied by the mole fraction to obtain a partial pressure for the 
carbon dioxide p^ 0 and water vapor p^ a . The mean beam length was estimated 
with the use of the table 4-2 in reference 14. The hot gases were assumed to have a con- 
stant density and to form a concentric cylinder about the plug. The factor by which the 
characterizing dimension (ref. 14) is multiplied to obtain was estimated by consider- 
ing two radiating gas configurations: (1) a rectangular parallelepiped 1 : °° : °° (infinite 
parallel planes) running parallel to the plug surface radiating to a spot on the plug sur- 
face (fig. 22(b)) and, (2) an infinite cylinder of half circular cross section running along 
the plug surface and radiating to a spot on the plug surface (fig. 22(c)). It can be seen 
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that the actual hot -gas configuration is less than the infinite parallel plane configuration 
but greater than the infinite cylinder of half cross-section. The factor was therefore 
estimated accordingly to be 1. 3. 

The characterizing dimension was determined by finding the perpendicular distance 
from the plug surface to the outer surface of the hot -gas jet at a large number of stations 
along the plug and averaging these values. For the hot -gas jet configuration at takeoff 
conditions, the characterizing dimension was calculated to be 0.427 meter (1. 4 ft). The 
mean beam length then became 0. 544 meter (1. 82 ft). The partial pressure of the 
COg, or of the water vapor, whichever is applicable, was multiplied by the mean beam 
length to give p^L^ or Using these factors and the average hot -gas static tem- 

perature, and e„„ at each station were determined from the charts of reference 14. 

The total -pressure corrections C CQ and C wa were taken to be equal to 1. 0 

(ref. 22). The corrections for the interaction of COg and HgO (Ae) were about 10 percent 

and were omitted. Since the correction reduces the radiation, this omission makes the 

wall-temperature prediction more conservative. The emissivity of the gas e p then 

equals + € (see eq. (C6)). Since an average hot-gas static temperature and aver- 
co wa, 

age hot-gas partial pressure of COg and HgO were used, e p is the same for any plug 
cross section. On the other hand, a^ w is different at each station since it depends on 
T w . Values of a CQ and a were obtained from charts (ref. 14) as functions of T w 
and Pcoh>( T w^p) or Pwa L^(T w /tp ). The values were corrected by multiplying by 

(V T w)°- 54 °r (tp/Tw) 0 ' 45 as applicable (eq. (C7)) to obtain a CQ and a^, respec- 
tively. Again, C CQ and were assumed to equal to 1. 0, and Aa was found to be 

about 4 percent and was neglected. In effect, this omission acts to increase the radiation 
and offset the effect of omitting Ae. 

The plug surface was assumed to be gray. A factor, (e + l)/2, to correct for this, 
was included in equation (C 2), which was developed for black bodies. The wall emis- 
sivity e was estimated to be 0. 65 based on data in reference 23. The effect of increas- 

vv 

ing the emissivity e w to 0.9 was also investigated (see fig. 14). 


Wall-to-Coolant Forced -Convection Heat-Transfer Rate ^ 

The convective heat transfer at the wall (eq. (C3)) was based on the assumption that 
the potential for heat transfer was the difference between the wall temperature and the 
adiabatic wall temperature T^ of the cooling film in a manner similar to that proposed 
in reference 24. 

The adiabatic wall temperature T^ at each station was determined with the use of 
the adiabatic wall film-cooling correlation proposed in reference 9 and modified to use 
an average value of the heat -transfer coefficient parameter (7rhD pI )^ from the pbint of 
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injection to the point of interest (as suggested in ref. 25). An experimental investiga- 
tion on the application of the preceding correlation to the surface of a plug nozzle is 
reported in reference 6. The correlation is repeated here: 


where 


In 


T 

rec, p 
T 

V rec, p 



fr h f D pl )o 

"c c p,sl 


0. 04 




h 


f 


= 0 


0265 


D 


L ReJ' 8 Pr 
Pi 


0.3 

f 



+ 0.4 tan 


-1 





1.0 


(C 8) 


(C9) 


(CIO) 


For this study, was always greater than V gl . 

The primary flow field properties were determined as indicated in the section on in- 
ternal flow properties. The local hot -gas recovery temperature Trp was determined 
assuming turbulent flow (ref. 20): 



(Cll) 


All coolant slot parameters (t gl , w g , Cp g j, q^, and V g j) were evaluated at slot exit 
conditions. 

A local heat -transfer coefficient h f x (eq. (C9)) was determined using air fluid 
properties evaluated at a temperature t f equal to an average of the local primary static 
temperature and the coolant slot inlet static temperature. The Reynolds number was 
based on 


Re f = PfV P D P 1 (C12) 

Mf 

where Vp is the local primary velocity. 

The average value of the heat-transfer parameter (^hf-Dpi)^ at a station was deter- 
mined by multiplying 77 times the local values of h f x and D pl and then averaging this 
parameter with those calculated at all local stations from the slot to the station of 
interest. 
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The local heat -transfer coefficient hg was calculated using the standard flat -plate 
correlation for turbulent flow (ref. 20) 


h 2 = 0.0296 -Re°' ^Pr 1 / 3 
X 


(C13) 


where fluid properties were evaluated at an average of the adiabatic wall temperature 
T w and the wall temperature T . The coolant flow boundary layer was assumed to be 
turbulent and to originate at the coolant slot (fig. 6(b)). The local velocity, for the cool- 
ant stream Reynolds number, was determined as indicated in the section on internal flow 
properties. 


Wall -to-Ambient Radiation (q^) 

The plug wall radiation is given by equation (C4) and depends on a configuration fac- 
tor and t^. The configuration factor Fg was shown in reference 26 to equal e w for 
gray bodies. Here the constant values of e = 0. 65 and 0. 9 were assumed over the plug 
length. The ambient radiation temperature t was assumed to be 311 K (560° R). 
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